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1 Introduction 

Many problems in science and engineering require solving nonlinear equation 

fi^) = 0, (1) 

[1-13]. One of the best known and probably the most used method for solving the preceding 
equation is the Newton's method. The classical Newton method is given as follows (NM) 

Xn+i = Xn - ^yp^, 71 = 0,1,2,3,..., and |/'(x„)| 7^ 0. (2) 

The Newton's method converges quadratically [1-13]. There exists numerous modifications of 
the Newton's method which improve the convergence rate (see [1-21] and references therein). 
This work presents a new quadrature based scheme for constructing optimal iterative methods 
of various convergence orders. According to the Kung-Traub conjecture an optimal iterative 
method based upon n + 1 evaluations could achieve a convergence order of 2". Through the 
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scheme, we construct optimal fourth order and eighth order iterative methods. Fourth order 
method requests three function evaluations while the eighth order method requests four function 
evaluations during each iterative step. The next section presents our contribution. 



2 Quadrature based scheme for constructing iterative meth- 
ods 

Our motivation is to develop a scheme for constructing optimal iterative methods. To construct 
higher order method from the Newton's method we use the following generalization of the 
Traub's theorem (see Theorem 2.4] and 13, Theorem 3.1]). 

Theorem 1. Let gi{x), g2{x),. . .,gs{x) he iterative functions with orders ri, ri,...,Ts, respec- 
tively. Then the composite iterative functions 

gix) = 9ii92i---igzix))---)) 

define the iterative method of the orders rir2r3 • • • r^. 

From the preceding theorem and the Newton method ([U, we consider the fourth order 
modified double Newton method 

_ f{x,i) 



Since the convergence order of the double Newton method is four and it requires four evaluations 
during each step. Therefore, according to the Kung and Traub conjecture, for the double 
Newton method to be optimal it must require only three function evaluations. By the Newton's 
theorem the derivative in the second step of the dobule Newton method can be expressed as 

f'{yn) = nx^)+ r r{t)dt, (4) 

J Xn 

let us approximate the integral in the preceding equation as follows 



fit) dt = UJi f{Xn) + i02 fiVn) + OJ3 /'(x„). (5) 

To determine the real constants wi, W2 and W3 in the preceding equation, we consider the 
equation is valid for the three functions: /(t) = constant, f{t) = t and t{t) = t^. Which yields 
the equations 

UJ1+LO2 = 0, 

Ll>i Xn +UJiyn+UJ3 = 0, (6) 

+ a;2 J/^ + ^3 2 a;„ =2(2/„-a;„). 

Solving the preceding equations and substituting the values in the equations ([4]) and (O, we 
obtain 

(7) 



/'(y„) = 2 (" /fa") 
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Combining the double Newton method and preceding approximation for the derivative, we 
propose the method (M-4) 

.f{Xn) 



o / JKVn) - f[Xn) 



fiVn) (8) 

f'iXn) 



\ l/n Xyi 

Since the method (M-4) is fourth order convergent and it requests only three evaluations. Thus 
according to the Kung-Traub conjecture it is an optimal method. We prove the fourth order 
convergent behavior of the iterative method ^ through the following theorem. 

Theorem 2. Let ^ he a simple zero of a sufficiently differentiable function / : D C R > R m 
an open interval D. If xq is sufficiently close to 7, the convergence order of the method (|8]) is 
4 and the error equation for the method is given as 

Here, e„ = — 7, Cm = /™(7)/™! with m > 1. 

Proof. The Taylor's expansion of f{x) and f'{xn) around the solution 7 is given as 

f{xn) = cie„ + C2el + C3el + c^e\ + O (e^ , (9) 
/'(x„) = Cl + 2 C2e„ + 3 c^el + 4 046^ + O (e^) . (10) 

Here, we have accounted for 7(7) — 0. Dividing the equations (|9]) and ()10|) we obtain 



f{Xr,) _ _ C2 2 o C3C1 - C2^ 3 3C4C1^ - 7C3C2C1 +4C2^ 4 , C 5\ 



From the first step of the method and the equations © and pUj) . we obtain 

= 7 + ^e„^ + 2 ^^^l^e„3 + - "™ + + O (e„^) . (12) 

Cl Cl^ Cl"^ 

By the Taylor's expansion of f{yn) around x„ and using the fist step of the method ([S]), we get 

/(...) ^ fM + /'(.,.) (-^) + (-||r^)' + ■ ■ ^ . (13) 

the successive derivatives of f{xn) are obtained by differentiating (ITOl) repeatedly. Substituting 
these derivatives and using the equations (jlip into the former equation 

/./ ^ 2 , ^ C3C1 - C2^ 3 3c4Ci^ - 7C2C3C1 +5c2^ 4 , 5\ n a\ 

/ 2/«) = C2e„ + 2 + ^ e„ + O (e^) . (14) 

Cl Ci^ 

Finally substituting from the equations ([9]) , (ITOl) and ([14]) into the second step of the contributed 
method (l8|), we obtain the error equation for the method 

e„+i = 3 e„ + O (ey . (15) 

Cl 

Therefore the contributed method ([8|) is fourth order convergent. This completes our proof. □ 
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To construct optimal eighth order optimal method, we consider the method 

f{Xn) 



Vn 



Vn 



fiVn) 



2 f fiVn) - fjXn 
fjZn) 



f'iXn) 



(16) 



f'(Zn)^f'{Xn)+ f"{t)dt, 

furthermore let us approximate the integral as follows 

J"{t) At = i/i /(X„) + 1/2 liVn) + f{Zn) + f'{x„), 



Since the order of the method ^ is four and order of the method ([2]) is two. Therefore by the 
theorem ([T]) convergence order of the method which is a combination of the methods 
and is eighth. The method ([TS]) require five function evaluations therefore, according to 
the Kung-Traub conjecture, it is not an optimal method. To develop an optimal method let us 
again express the first derivative by the Newton's theorem 

(17) 



(18) 



to determine the real constants, i^i, V2, 1^3 and in the preceding equation, we consider the 
equation is valid for the four functions: /(t) = constant, /(t) = t, t{t) — and /(t) — t^. And, 
we obtain the four equations 

1^1 + 1^2 + 1^3 = 0, 

Xn + 1/2 Vn + VzZn + Vi = 0, 

Vl xi + l'2Vn + ^3^1^ + 2 ViXn =2 (z„ - Xn), 

vixl + V2vl + v^zl + ivixl =2,{zl-xl). 
from the preceding equation and the equations (flT)) . ([18]), we get 

1 



f'iZn) - 



{-Vn+Xn) {- Zn + Vn) {- Zn + Xn) 



{-Zn + Vn) {-Zn + Xn) {-Vn + Xn) f'{Xn) 



(Xn-Vn) {2Vn - SZn + Xn) f (Zn) + {Xn - Zn) / (j/n) - (j/ri - Z„) (3 a;„ - 2 y„ - Z„) / (a;„) 



Combining the eighth order method (|16p and the preceding equation, we propose the following 
optimal eighth order iterative method (M-8) 

f{Xn) 



Vn 



f'iXnV 



Vn 



fiVn) 



fiVn) - f{Xn) 



l/n X n 



f'iXn) 



fiZn) 



{-Vn+Xn) {-Zn+ Vn){-Zn+ Xn) 



(19) 



{-Zn + Vn) {-Zn + Xn) 



{-Vn + Xn) f'{Xn) - {Xn - Vn) (2 Vn " 3 Z„ + Xn) f {Zn) 
+ {Xn - Zn)^ f {Vn) - {Vn " Zn)'^ (3 Xn - 2 y„ - Z„) / {Xn) 
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Since the method (M-8) is eighth order convergent and it requests only four evaluations during 
each iteration. Thus according to the Kung-Traub conjecture it is an optimal method. We 
prove the eighth order convergent disposition of the iterative method ^ through the following 
theorem. 

Theorem 3. Let j be a simple zero of a sufficiently differentiable function / : D C R i— > R m 
an open interval D. If Xq is sufficiently close to 7, the convergence order of the method (jl9p is 
8. The error equation for the method (I19p is given as 

C2^ (C3C1^C4 - C4Ci2c2^ - Ca^Ci^Ca +2C3C1C2^ - C2^) g ^, ,9 

en+i = ^ ^ -e„ + O (e„) . 

Proof. Substituting from the equations ([9]), (fTO|). (fTl]) . (fT4|) into the second step of the con- 
tributed method (fT9l) yields 



= 7 - ^^^^IZ^^l _ 2 ^2C4C1^+C3^C,^-4C3C,C2^ + 2C2^ ^^, _^ ^ ^^^^ 
Ci Ci 

Here, we have used the first step of the method ([T9|. To find a Taylor expansion f{zn), we 
consider the Taylor's series of f{x) around yn 

f{Zn) = fiVn) + f'iVn) {Zn ^ Vn) + ^^^{Zn " Vn)' + ' ' ' , (21) 

substituting from equation ([Ti)) and using the second step of the contributed method we 
obtain 

/(.„) = - ^"''^ "' en' - 2 ''''''' + - '^^^^^^^ + ^'^^' e..^ + O . (22) 

Here, the higher order derivatives of f{x) at the point y„ are obtained by differentiating the 
equation ([Tl)) with respect e„. Finally, to obtain the error equation for the method p^ . substi- 
tuting from the equations dH), ([10]), (HH), ([121), dSHl) ([221) into the third step of the contributed 
method ()19p yields the error equation 

C2^ (c3Cl3c4 - C4Ci2c2^ - C32ci2c2 -H2C3C1C2^ - C2^) o ^ , ,9 

en+1 = ^^^^^ el + O (e„) , (23) 

which shows that the convergence order of the contributed method (|TO)) is 8. This completes 
our proof. □ 



3 Numerical examples 

Let us review some well known methods for numerical comparison. Based upon the well known 
King's method [12] and the Newton's method recently Li et al. constructed a three step 
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and sixteenth order iterative method (LMM) 



Vn 



Vn 



2f{Xn)-f{yn) fiVn) 



Xn+1 



f 



f{Zn) 



(24) 



fiZn) 



[l6| . Based upon the Jarratt's method fit], recently Ren et al. [iTlfli 
convergent iterative family consisting of three steps and two parameters (RWB) 

2 f{Xn) 



formulated a sixth order 



Vn = 




Zfi — 









3/'(a;„)' 

3f(yn) + f f{Xn) 

6/'fy„)-2/'(a:„) fixr^y 

(2a-6)/'(a;„) + 6/'(y„) + c/(x„) 



(25) 



f{Zn) 



(-a - h)r{xn) + (3a + 6)/'(2/„) + cf{x,,) f'{xn) ' 

where a, &, c G M and 07^0. Wang et al. [18] also developed a sixth order convergent iterative 
family, based upon the well known Jarratt's method, for solving non-linear equations. Their 
methods consist of three steps and two parameters (WKL) 

2 f{Xn) 



y-n 



Xn+1 = 



3/'(a^„)' 

3f'iyn) + f'{Xn) f{Xn) 



6/'(y„) - 2/'(x„) f'ixnY 

(5a + 3/3)/'(x„) - (3a + /3)/'(2/„) /(z„) 



(26) 



2a/'(a:„)+2/3/'(y„) 



where a, /3 e M and a + /^t^O. Earlier, Neta (20[ has developed a sixth order convergent family 
of methods consisting of three steps and one paramete (NETA) 

nxnV 

fjXn) + afjyn) fiVn) 

/(x„) + (a - 2)/(2/„) /'(x„)' 

fjXn) - fiVn) f{Zn) 
f{Xn) - 3/(y„) f'(Xn)' 

We may notice that, in the preceding method, with the choice a — —I the correcting factor in 
the last two steps is the same. Chun and Ham [2l[ also developed a sixth order modification 
of the Ostrowski's method. Their family of methods consist of three-steps (CH) 

fi^n) 



Vn = 




Zfi — 


Vn 


Xn+1 ^ 





yn = 


Xn 


Zfi — 


yn 


Xn+1 = 





fjXn) fjyn) 

/(a;„) - 2/(y„) /'(x„)^ 



(27) 
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where u„ — f{yn)/f{xn) and'H(t) represents a real valued function satisfying ?^(0) = 1, H'{0) = 
2. In the case 

the third substep is similar to the method developed by Sharma and Guha [l^. The classical 
Chebyshev method is expressed as (CM) 



f'iXn) 



2 /'(x„)2 



|l5j and the classical Halley method is expressed as (HM) 



2r{Xn?-f"{Xn)f{Xn) 



(29) 



(30) 



|15| . The convergence order ^ of an iterative method is defined as 



hm I—— =C7^0, 



and furthermore this leads to the following approximation of the computational order of con- 
vergence (COG) 

_ ln|(x„+i - 7)/(x„ - 7)1 
In |(a;„ - 7)/(a;„_i - 7)] ' 

For convergence it is required: \xn+i — Xn\ < e and |/(a::n)| < £■ Here, e — lO^"^^*^. We test the 
methods for the following functions 



/i(^ 

f2{x 



= +4a;2 - 10, 

— X exp(x^) — sin^(x) + 3 cos(x) + 5, 
= sin^(x) - + 1, 

= tan^"'^(a;) 

— x'^ + sin(7r/a;^) — 5, 

= e(-"'+"+2) - 1, 



7W1.365. 
7w - 1.207. 
7« ± 1.404. 
7=0. 
7=\/2. 
7=2.0. 



Computational results are reported in the Tableland the Table [2] The Table [T] presents (num- 
ber of functional evaluations, COC during the second last iterative step) for various methods. 
While the Table [2] reports \xn+i ~ Xn\ for the method M-8. Free parameters are randomly 
selected as: for the method RWB a = 6 = c = l, in the method by Chun et al. (CH) (3 = 1, 
in the method WKL a = /3 = 1, in the method NETA a = 10. 

An optimal iterative method for solving nonlinear equations must require least number of 
function evaluations. In the Table[Tl methods which require least number of function evaluations 
are marked in bold. We acknowledge, through the Table [1] that the contributed methods in 
this article are showing better performance to the existing methods in the literature. 
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Xq 


HM 


CM 


LMM 


NM 


RWB 


NETA 


CH 


WKL 


M-4 


M-8 


/l(^) 


1.2 


(27, 3) 


(27,3) 


(24, 16) 


(20, 2) 


(20,3) 


(20,6) 


(20,6) 


(20,3) 


(18,4) 


(16,8) 


M^) 


-1.0 


(24,3) 


(24,3) 


(24,15.5) 


(22, 2) 


(20,3) 


(20,6) 


(20,6) 


(20,3) 


(18,4) 


(16,8) 


M^) 


1.5 


(21,3) 


(21,3) 


(18,15.8) 


(20, 2) 


(20,3) 


(20,6) 


(20,6) 


(20,3) 


(18,4) 


(16,8) 




0.5 


(21,3) 


(21,3) 


(24, 24) 


(18,2) 


(20,3) 


(16,7) 


(16,7) 


(20,3) 


(18,5) 


(16,11) 




1.3 


(24,3) 


(24,3) 


(24, 16) 


(20, 2) 


(20,3) 


(20,6) 


(20,6) 


(20,3) 


(18,4) 


(16,8) 




1.2 


(24,3) 


(27,3) 


(24, 16) 


(26, 2) 


(20,3) 


(20,6) 


(20,6) 


(20,6) 


(21,4) 


(20,8) 



Table 1: (number of functional evaluations, COC) for various iterative methods. 





f2{x) 


Mx) 


Mx) 




Mx) 


Mx) 




1.6 X 10"^ 


2.0 X 10"^ 


9.5 X 10"^ 


4.2 X 10" 


1 


1.1 X 10"^ 


7.9 X 10" 


1 


3.3 X 10"'' 


2.4 X 10"® 


4.0 X 10"^° 


1.7 X 10" 


5 


1.1 X 10~* 


8.6 X 10" 


4 


6.4 X 10"'^^ 


1.9 X 10"*^ 


6.3 X lO"'^'' 


3.1 X 10" 


55 


2.8 X lO"**^ 


2.9 X 10" 


25 


1.1 X 10"'^*'^ 


2.8 X lO"^^** 


2.5 X 10"''" 


1.9 X 10" 


601 


5.8 X 10"^^* 


5.8 X 10" 


197 












1.3 X 10" 


1570 



Table 2: Generated \xn+i ~ x„ \ with n > 1 by the method M-8. For initialization see the 
Table m 
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